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Abstract 

Using the vortex filament model and the Gross Pitaevskii nonlinear Schroedinger equation, we 
show that bundles of quantised vortex lines in helium II are structurally robust and can reconnect 
with each other maintaining their identity. We discuss vortex stretching in superfluid turbulence 
and show that, during the bundle reconnection process, Kelvin waves of large amplitude are gener- 
ated, in agreement with the finding that helicity is produced by nearly singular vortex interactions 
in classical Euler fiows. 

PACS numbers: 67.25.dk Vortices and turbulence, 
47.27.De Coherent structures, 
47.32.cf Vortex reconnection and rings 



Concentrated, tubular vortical regions are commonly observed in turbulent flows. These 
vortex structures are important in turbulence dynamics, for example they enhance mixing 
and diffusion. The role which they play in determining turbulence statistics (spectra, in- 
termittency) is the subject of much study [1]. If the axes of tubular vortex structures are 
interpreted as the skeleton of turbulence, then the knottedness of the axes characterises 
the topology, and vortex reconnectionsQ] are the critical events which change this topol- 
ogy. This idealised picture becomes reality if one moves from ordinary fluids to superfluids 
(liquid ^He and ^He, atomic Bose-Einstein condensates and neutron stars). The reason is 
that in superfluids quantum mechanics constrains any rotational motion to ultra thin vortex 
filaments whose core radius and circulation are fixed, and turbulence takes the form of a 
tangle of such discrete filaments {3]. Recent experimental, theoretical and numerical studies 
have revealed similarities between ordinary turbulence and superfiuid turbulence 4l], such 
as the same Kolmogorov -5/3 energy spectrum [5] and the same temporal decay law of the 
vorticityp]. 

This letter is concerned with vortex reconnections and vortex stretching in superfiuid 
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turbulence. The importance in helium II of vortex reconnections and its scaling laws[7|] 
was first appreciated by Schwarz[8], who developed the reconnecting vortex filament model; 
later, the existence of reconnections was confirmed by Koplik and Levine[!|] using the non- 
linear Schroedinger equation (NLSE) model. Vortex reconnections are associated with the 
dissipation of superfiuid kinetic energy in the limit of absolute zero lOj, either directly HI] or 
via a Kelvin wave cascade[l2|. Recently, individual vortex reconnections have been detected 
in experiments |14l] . As recognised by Procaccia and Sreenivasan[15], the occurrence of the 
-5/3 law in superfiuid turbulence is surprising if one notices that vortex stretching, usually 
recognised as an important mechanism to transfer energy across length scales, is absent 
because the radius of the superfiuid vortex core is fixed. A possible solution of the puzzle is 
that some filaments are organised in vortex bundles, which have indeed been noticed in the 
most recent numerical simulations of superfiuid turbulence IGj. Clearly, stretching occurs 
if the relative position of vortex strands within a bundle changes during the evolution. 

In a first set of numerical experiment we use the model of Schwarzjs]. Let the space 
curve s(^, t) represent a vortex filament where ^ is arclength and t is time. The self-induced 
velocity v^j of the vortex at s is 
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where k is the quantum of circulation and the Biot-Savart integral is suitably 
desingularized 8'1 . The numerical method to evolve configurations of vortex filaments has 
been descHbed elsewhere gQ. Here ,t suffiees to .a. that the ma,„e.ts a.e disc.et.sed ,nto 
a large, variable number of points, A^, depending on the local curvature: points are removed 
in regions where filaments straighten and are added where the local radius of curvature 
becomes smaller. The time evolution is based on a 4*^ order Runge-Kutta scheme with 
variable time step At, which depends on the current minimum distance 5min along points, 
which determines the frequency of rotation of the fastest Kelvin wave perturbation. If two 
vortex strands become closer to each other than the local discretization along filaments, 
then, consistently with the orientation of the filaments, our numerical algorithm reconnects 
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the strands [17], provided that the total length is decreased [18|]. This condition [19|] mimics 
the small kinetic energy losses at vortex reconnections discussed in Ref. Jll|. We also checked 
that, with our algorithm, the reconnection condition used by Schwarz[8] is satisfied [2^. 

Using this model, we study the interaction of two bundles of a given number M of 
(initially) straight parallel vortex strands, set (initially) at 90 degrees to each other. Le D 
be the distance between the axes of the two bundles and A be the radius of each bundle. 
The calculation is performed in a cubic periodic box —B < x,y, z < B. A typical result for 
M = 7, A = 0.0417 cm, D = 3A and B = 1 cm is shown in Fig. [H The initial position of 
vortex strands within the same bundle is symmetric: in this case we place six vortices at the 
corner of a hexagon and one vortex in the middle. In the absence of the the second bundle, 
the vortex strands of the first bundle would rotate around each other in a time scale r of the 
order of 10 s because they are parallel. With two bundles, the interaction makes them to 
bend in the direction of each other, until, at time ti ^ 6.7 s, the first reconnection occurs. 
The successive evolution involves the reconnections of all strands at times ^2 ~ 16.4 s, 
^3 ^ 21.0 s, t4 ^ 23.6 s, ts ^ 27.4 s, ~ 52.0 s, and ^ 65.9 s, after which the two bundles 
separate from each other (Fig. [1], bottom right) and move away. 

Similar calculations with different parameters show that the last reconnection (at which 
the remaining vortex bridge breaks, and the bundles become free from each other) usually 
takes longer than the first few reconnections. Calculations also show that, after the first re- 
connection, the vortex strands develop Kelvin waves of large amplitude. During the process. 



the total vortex length L increases by about 40%, as shown in Fig. [2l Fig. [3] shows the aver- 
age inverse radius of curvature, < 1/R>, obtained by computing \(fs/d^'^\ at each point Sj 
(j = 1, ■ • ■ A^) and then averaging over all points. To cope with the increase in L and decrease 
of < i? >, the number of points (initially = 700) grows with time, up to = 4191, when 
we stop this particular calculation. Fig. H] shows normalised histograms (Probability Density 
Functions) of the average inverse curvature computed at different times: the formation of 
the Kelvin wave cascade is evident in the development of a tail of the distribution at larger 
and larger times. All these results were confirmed by further calculations with M = 5 and 
11 strands. 

To make sure that our result does not depend on the reconnection algorithm, we study 
the interaction of vortex bundles by solving the NLSE, also called the Gross-Pitaevksii (GP) 
equation: 

2^^ = -VV + |^|V-^- (2) 



The NLSE is a convenient model of superfluid helium II 2l|; for example, Koplik and Levine 
9(] used it to confirm Schwarz's insight that quantised vortices reconnect. Eq. [2]is written in 
terms of h/\/2mE (the coherence length) as unit of space, 2E/h as unit of time and mE/V 
as unit of density, where m is the mass of one boson, V the strength of the interaction 
between bosons, and E the energy per boson. The calculation is done using a 5^^ order 
Runge Kutta Fehlberg method (with typical time step of the order of 0.01) using 256^ grid 
points and reflective boundary conditions. 

In a second set of calculations, we thus solve the NLSE in a cubic box —128 < x,y, z < 
128. The initial condition consists of two bundles of radius A = 8 and distance D = 3A 
between the axes of the bundles. Each bundle contains M = 7 vortices (six at the corners 
of a hexagon and one in the middle). The time sequence, shown in Fig. [5l confirms the 
previous result that vortex bundles maintain their identity and reconnect with each other. 
Note the emission of small vortex rings in the last image. As in the previous calculation, 
large amplitude Kelvin waves are generated. Fig. E] shows that at first L increases with 
t, then, when the bundles have moved sufficiently away from each other, L saturates; the 
relative increase of L is about 30%, confirming the intensification of vorticity. The relation 
between energy and length is important; it is common in the literature to state that the 
energy per unit length of vortex hne is (p-jK^/ (47r)) In {b/a) where ps is the superfluid density 



and b is an upper cutoff (the radius of the container or the distance to the next vortex), 
but this relation is vahd only for a straight vortex of course. The NLSE model allows to 
determine the energy more precisely than the vortex filament model: in our calculation 
the total mass and the total energy are conserved within 1 part in 10^ and 1 part in 10^ 
respectively. Further calculations with M = 5 and 9 strands confirm our results. 

The above results refer to temperature T = 0. If T > the motion of a vortex filament 
is governed by Schwarz's equation 

ds 

— = v^i + as' X (v„ - Vsi) - a's' x [s' x (v„ - v^j)], (3) 

where s' = ds/d^ is the unit tangent at s, v„ is the normal fluid velocity and a and a' 
are known temperature dependent coefficients resulting from the mutual friction between 
the vortex lines and the normal fluid 22|. In a third set of calculations we find that, 
at T = 1.65 K (the typical temperature of many experiments), with v„ = , bundle 
reconnections are still possible, although with much less helical disturbances, as shown in 



Fig. [71 in agreement with calculations [17| which show that at T > the vortex tangle is 
much smoother than at T = 0. We do not repeat this calculation with the NLSE model 
because there is not yet a consensus on how to generalise the NLSE to finite temperatures. 



23|. 



although many approaches have been proposed 

We conclude that vortex bundles are structurally stable structures, in the sense that they 
survive a time longer than their characteristic time r, travel a distance larger than their 



size A, consistently with results for bundles of vortex rings 2J]. Remarkably, vortex bundles 
survive reconnections with other bundles without disintegrating, but rather amplifying their 
vortex length. Finally, the coiling of the vortex strands which is triggered by bundle recon- 
nections confirms results of Holm and Kerr 25|] about the generation of helicity in nearly 



singular vortex interactions of the Euler equation. 
C.F.B. is indebted to W.F. Vinen for discussions. 
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FIG. 1: (Colour online). Reconnection of two bundles of seven vortex strands each, a, top left: 
t = s. b, top right: t = 7.13 s. c, middle left: t = 23.58 s. c, middle left: t = 36.27 s. d, bottom 
left: t = 61.49 s. e, bottom right: t = 80.35 s. 
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FIG. 2: (Colour online). Total vortex length L vs time t corresponding to Fig. [TJ 
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FIG. 3: (Colour online). Average inverse radius of curvature < > vs time t corresponding to 
Fig.ffl 
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FIG. 4: (Colour online). PDF of the inverse radius of curvature at different times, corresponding 
to Fig. [H From left to right, the times are: (solid red): t = 7.13 s, (solid green) t = 23.58 s, 
(dashed blue) t = 36.27 s, (solid purple) t = 61.49 s, (sohd hght blue) t = 80.35 s. 
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FIG. 5: Reconnection of two bundles of seven vortices each computed by solving the dimensionless 
NLSE. The figures show density isosurfaces at the level 0.25 (of the unit bulk density away from 
vortices) at different times t. The numerical resolution of the vortex core is such that at 0.25 
density level there are about 3 grid points within a vortex core, and at 0.90 density level there are 
about 15 grid points. Top right: t = 110; Middle left: t = 240; Middle right: t = 320; Bottom left: 
t = 450; Bottom right: t = 800. 
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FIG. 6: Dimensionless vortex length L vs dimensionless time t corresponding to Fig. [5l 
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FIG. 7: Comparison between bundles reconnections at T = K (left, corresponding to a = a' = 0) 
and at T = 1.65 K (right, corresponding to a = 0.111 and ol = 0.01437). The initial condition is 
the same for both calculations. The times (left: t = 86.9 s, right: t = 58.1 s) are chosen so that in 
each case the bundle reconnection has proceed to the point that only two vortex strands are still 
part of the initial bundle. 
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